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ABSTRACT: The nucleoli effective mass, binding energy and neutron star configura- 
tion in the Relativistic Mean-Field Theory (RMF) is considered. The calculation are 
motivated by the construction of the equation of state for a neutron star in the RMF. 
The equation of state for the parameters set TM1 is calculated using the Feynman 
- Bogolubov variational method for temperature different from zero. The structure 



of the neutron star is presented. The maximal stable configuration are obtained for: 
M max = 1.91 M , R = 12.84 km. 
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Introduction 



The physics of compact objects like neutron stars offers an intriguing interplay be- 
tween nuclear processes and astrophysical observables. Neutron stars exhibit condi- 
tions far from those encountered on earth. The determination of an equation of state 
(EoS) for dense matter is essential to calculations of neutron star properties. 
This paper presents a basic model of neutron star (l| matter including interactions 
among nucleons in the relativistic mean field approximation @@[|]. Especially the 
Walecka model (QHD) and its nonlinear extensions have been quite successful and 
widely used for the description of hadronic matter and finite nuclei. Increasing in- 
terest in neutron matter at finite temperature has been observed recently in relation 
to the problems of hot neutron stars and of protoneutron stars and their evolutions 
in particular. Theories concerning protoneutron stars are being discussed in works 
by Prakash et. al. ||. Recently a detail calculation has been done with different 
models to study the properties of neutron stars ||. Glendenning @ has studied the 
properties of neutron star in the framework of nuclear relativistic field theory. In our 
calculations, we used the TM1 [[<]] parameter set, which has a capability to reproduce 
the known results of finite nuclei as well as of normal nuclear matter. The TM1 
model possess a distinctively stiffer EOS. 

The Relativistic Mean Field Theory 

The fields of the model RMF for a, u and p-mesons are denoted as <p, lu^, p M . The 
Lagrange density function for this model has the following form 

c = \d^d^ - \r%r^ - \f, u f^ + \mI ^ + \mI P y» (i) 

-U{<p) + ^c 3 (^) 2 + i^D M V - i>{M - g a <p)i/> + 

2 2 

i £ TTffdpLf - J2 9f(L f He Rf + h.c.) (2) 
/=i /=i 

where 

R% = d,p a u -duPl + ge abcP l P c u (3) 

Fp, = d^ujy - dyU)^ (4) 

= d» + -^9 P Py a + i9ui^i* (5) 

The potential is given by 

U((p) = l -m 2 s ip 2 - ^ 2 <p 3 - ^ 3 <p 4 = ^m 2 <p 2 + ^K<p 3 + ^A<p 4 (6) 
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Parameter 


L2@ 


NBL0 


LN1 12 


TM10 


M 


938 MeV 


938 MeV 


938 MeV 


938 MeV 


M w 


786 MeV 


786 MeV 


795.359 MeV 


783 MeV 


M p 


770 MeV 


770 MeV 


763 MeV 


770 MeV 


m s 


500 MeV 


510 MeV 


492 MeV 


511.198 MeV 


92 = k/2 





2.03 /m- 1 


12.172 fm- 1 


7.2325 /m- 1 


93 = A/6 





1.666 


-36.259 


0.6183 


9s 


10.0773 


9.6959 


10.138 


10.0289 


9u> 


13.8655 


12.5889 


13.285 


12.6139 


9 P 


8.48784 


8.544 


4.6322 


4.6322 


c 3 











71.3075 



Table 1: Parameters Set for the Lagrangian (||). 



The fermion fields are composed of protons, neutrons and electrons, muons and 
neutrinos 



Li 



e Rf = V.R; V>r) 



(7) 



M is the nucleon mass and m 8 , M u , M p are masses assigned to the mesons fields, g, g' 
and g s are coupling constants. The Lagrangian function includes also the nonlinear 
term |C 3 (o; (U a; M ) 2 which affects remarkably the form of the equation of state, the 
Higgs field H takes the form of 



takes here the residual form. The parameters used in NBL model 
MeV with k = 800 MeV and A = 10. The Euler equation for $ A 
fields are 

Dtp = m 2 s (p + g 2 (p 2 + g 3 ip 3 - g s ipip 



(8) 

are m s = 500 
{ip, p M ,V4 

(9) 



where 

is the baryon current 
and 



-D p R^ a = M 2 pP ^-g p J» 



(10) 

(11) 
(12) 

(13) 
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Figure 1: The potential U(ip) in the relativistic mean field theory for L2 0, NBL||, LN1 
fi~2||, TMiP] phenomenological parameters. 



is the isospin current. In the system we have conservation of baryon charge 

Qb = J d 3 x.J°, 

and the isospin charge 

Q 3 = J d 3 xJ° 3 . 

The last is the Dirac equation 

irfB^ — (M — g s <p)i> = 0. (14) 
The physical system is totally defined by the thermodynamic potential || 

n = -kTlnTr{e~ p{H -^ QB -^ Qa) ) (15) 
where H is the Hamiltonian of the physical system 

H = Y. I d 3 x{d <S> A ir£ - £} (16) 

A J 

and 7r A = d(^$ A ) is a momentum connected to the field The fields &a = {cp, lo^, 
p^,ip} denote all fields in the system. The average charges 

on ^ on 

■jt = ~ <Qb>, 7 — = -<Q 3 > 17 
a\i a/i 3 

can be obtained from the thermodynamic potential, of course they should be con- 
served. In this paper we shall use the effective potential approach build using the 
Bogolubov inequality |§ 

Q < fij = tt (m Eh m F )+ < H - H > (18) 
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fio is the thermodynamic potential of the trial system as effectively free quasiparticle 
system described by the Lagrange function 



C (m B , m F ) = \d»Tpd»Xp - \m\Tf? - -^G^ - \k., u F^ 

" 1 
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+ X -Mlw^ + l -M 2 p -p a ^ + H^D, - m F )^ (19) 



Similar to the general case 

G~l„ = d^-p a v - d„pl 

and 

Fu)^ v = d^uj v — <9„u7 M . 

We decompose the <&a held into two components, the effectively free quasiparticle 
field $,4 and the classical boson condensate £4 

$A = §A + U (20) 

In the case of the RMF model we have 

(f = Tp + a (21) 

uj^ = uJ^ + «V, = S^ow (22) 

P a »=r + r% r; = 5 a '%, r (23) 

The £a = w, r} field will be treated as the variational parameters in the effec- 
tive potential. Also the boson and fermion mass me, will be treated as as the 
variational parameters. The covariant derivative for the trial system is 

Df, = + ^ig p r a ^a a + igjw^ (24) 

This introduce the homogenous fermion interaction with boson condensate w^, r®. 
The fermion quasiparticle will obey the Dirac equation 

{ij' i D ll -m F )^ = (25) 

The constant condensate w, r simply shift the chemical potential from pi = p® (when 
w = r = 0) to 

V>n = + ^9 p t - g u w (26) 
where p n = p + \p$ and p p — p — |/x 3 . 

Neutrons, protons and electrons are in /^-equilibrium which can be described as a 
relation among their chemical potentials 

p p + jJ, e = jJ, n (28) 
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where /jl p , /i n and \i e stand for proton, neutron and electron chemical potentials 
respectively. If the electron Fermi energy is high enough (greater then the muon 
mass) in the neutron star matter muons start to appear as a result of the following 
reaction 

e" -> \T + v e + (29) 

The chemical equilibrium between muons and electrons can be described by the 
condition 

(J>[i = He (30) 

When neutrinos are trapped inside the protoneutron star also the neutrino chemical 
potential should be included 

H P + n e = fi n + n Ve . (31) 

The density of the thermodynamic potential f\ = Q±/V is equal to 

f 1 {m B ,m F ,a,w,r)= (32) 
fo(m B , m F ) + \ < Tp 2 > (m 2 s - m%) + |A < Tp 2 > 2 + 
\k <Tp 2 > a + \{m 2 s + \\<Tp 2 > )o- 2 + 
ly/to- 3 + ^Act 4 + < ipip > (go - m F ) + •-, 



fo = f B + f F (33) 

where fs is the boson free energy and fp the fermion one. For boson field the free 
energy is 

/B = ^/ rfVn(1 - e_MP)) (34) 



with oj(p) = y p 2 + Tn B . For fermions we have 4 degree of freedom, 2 for spin 1/2 
and 2 for particle - antiparticle distinguishing. For one fermion field the free energy 
is equal to 

fp = ~ E / 'd 3 p{ln(l + e-^-^) + ln(l + e-^ + ^ (35) 



i={n,p} 



now with e(p) = \Jp 2 + rn 2 F . Variation 



dm 2 B dmp 



with respect to the trial system L gives 

,2 _2 i o i 0„ ( _2 i ^ ^2 



m£ = mf + 2^ 2 <T + 3^ 3 (<r + < V? >o), (36) 
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Figure 2: Meson mass m 2 B /M 2 dependence on the Fermi momentum x = kp/M for 
temperature T = 50MeV. 

m F = MS = M — g s a. (37) 

In the local equilibrium inside the star the free energy reaches the minimum at a. 

The same result may be achieved calculating the averages of the equation of 
motions (H) for the effective system C . In the mean field approximation the meson 
field operators are replaced by their expectation values. We also consider the isotropic 
system at rest. For the scalar field the equation (||) have is follows 

1 11 — 1 

(m 2 s + -\ <Tp 2 > )a + -ko 2 + —Act 3 = g s < ijnf) > --K < > (38) 
I I o I 

If we shall notice that 



dfB 1 _ 2 ^ 

dm 2 , = 2 < V >0 ' 



then it easy to obtain 



< V 2 >o~- 



2tt 2 



dpp 2 



^Jp 2 + m 2 B {exp(Puijp)) 



This result means that equation ( j36l ) is highly nonlinear with respect to the meson 
mass ttib (Fig. ||). As < ipip > (calculated with respect to £ system) depends 
on the effective nucleon mass mp (or a ) the equation (|38|) is highly nonlinear also 
with respect to a. In most of papers [|1(J temperature dependence of the boson fields 
is neglected. Contrary to the fermion case the effective mass of boson m# growing 
with Fermi momentum. In the result the bosons temperature contributions may be 
neglected. 

Calculation similar to the boson case, base on the relation 

amp 
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Figure 3: The effective neutron mass mp = 5(xp)M as function of the neutron Fermi 
momentum x = kp ( /m _1 ). 



gives 



i={n,p} 



—r , x , m F f°° p 2 dp , 1 . . 

<w>o= .-_?^i» ^r^r { -p(/3fo-A,))+i + (39) 

1 ■}. 



exp(/3(e p + ft)) + 1 

The quantum average < ipijj >o depends on the neutron and proton chemical po- 
tentials fjLp, /i n . In the result the effective nucleon effective mass m F also will be 
dependent he neutron and proton chemical potentials. In the result the solution a of 
the equation (P) also will be dependent on £. The same situation will consider other 
fields. This model is the simple example of the relativistic mean filed theory 0. 

The effective mass mp (or 5 = mp/M ) dependence on the dimensionless Fermi 
momentum xp is presented on the Fig|| The binding energy 

E = e(x F ,T)/Q B -M (40) 

for nucleon symmetric phase is presented on Fig. 2 (see Table 0). We see that 
binding energy strongly depends on temperature and above T > 15 MeV is positive. 
The temperature dependence of the binding energy in presented on the Fig |j. To 
calculate the properties of the neutron star we need the energy-momentum tensor. 
In case of the fermions field it is more convenient to use the reper filed defined as 
follows g^ u = e a e b u r]ab where T) a b is the flat Minkowski space-time matrix. The general 
definition 

„ dLn „ OLp 

T >» = 2 dg^ + <d£ ~ 9 ^ L (41) 
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Parameter 


L2§ 


NBL0 


LN1|2| 


TMlflTJ 


E Q MeV 


-15.74 


-16.59 


-16.42 


-16.26 


kF,o fm' 1 


1.30 


1.31 


1.309 


1.29 


po fm~ 3 


0.149 


0.1524 


0.1517 


0.1452 


5q = rapjM 


0.54 


0.60 


0.573 


0.634 


K MeV 


545.93 


445.73 


215.19 


281.53 



Table 2: The nucleon symmetric phase properties at the saturation point. 



E/Q-M ( MeV ) 



T=25 MeV 



T=0 MeV 



TM1 




M fm ) 



0,5 



2,0 



Figure 4: The temperature dependence of the binding energy for nucleon symmetric TM1 
phase. 



allows us to calculate the density of energy and pressure. The total Lagrange function 
L = L B + L F is divided into boson and fermion part. To calculate the density of 
energy and pressure we shall average the energy-momentum tensor T^ v with respect 
to the quasi equilibrium configuration defined by the trial system L . We define the 
density of energy and pressure by the energy - momentum tensor 

<T^>=(P + e)u^u u - P 9lw (42) 

where is a unite vector (u^u^ = 1). So, the calculations give 

e(x F , T) = pc 2 = - l -Mlw 2 - ~M p V + g u Q B w + \g P Q 3 r - ^c 3 w 4 + U{a) + e F (43) 



1 1 1 

' - r2„..2 , 1 A/ r2 2 , 1 „ 4 



where 



P(x F , T) = -Miw 1 + -My + -c 3 w* - U(a) + P F 



£f = e x(xp,T) 

Pp = P 0( j)(xF,T) 



(44) 

(45) 
(46) 
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L2 NL1 NBL 



TM1 




100 200 300 400 500 600 700 800 900 1000 



Figure 5: The equation of state for different parameters sets (L2, NBL,TM1). 

The fact that neutron mass depends on fermion concentration ( or neutron chemical 
potential //) now must be included into x( x f,T) and <j)(x F ,T), 



1 f°° i 

x(x F ,T) = ^ dzz 2 ^z 2 + 5 2 (x F ){ 



1 



+ 



exp((^j6 2 (x F ) + z 2 - n')/r) + 1 
1 



exp((J5 2 (x F ) + z 2 + fi')/r) + 1 



(47) 



{x F ,T) 



z A dz 



-{- 



1 



3tt 2 7o ^z 2 + 8 2 {x f ) exp(( v /(5 2 (x F ) + z 2 - h')/t) + 1 

1 



exp((^ 2 (a; F ) + z 2 + /x')/r) + 1 



} 



where r = (k B T)/M, 



and 



H 1 = nfM = J5 2 {x F )+x 2 F 



x F = k F /M 



(48) 



(49) 



(50) 



Similar to paper |T2| we have introduced (|49|j50|) the dimensionless "Fermi" momen- 
tum even at finite temperature which exactly corresponds to the Fermi momentum 
at zero temperature. Both e F and P F depend on the neutron chemical potential /x or 
Fermi momentum x F . This parametric dependence on /i (or x F ) defines the equation 
of state. The various equations of state for different parameters sets is presented on 
fig|| The equation of state for the parameters set TM1 for temperature T=10 MeV 
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T=20 MeV 




P'Po 



0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 

Figure 6: The equation of state for the parameters set TM1 for T=0 (solid line) and 
T=20, 50 MeV temperature. 

is presented on fig|| It is interesting to notice that even for Qb = due to the 
presence of the thermal exited particle antiparticle pairs there is finite energy and 
pressure density. It is interesting to notice that the meson field <p as the scalar field 
contributes to the pressure with negative sign while the vector meson fields (u, p) 
with the positive one. 



The neutron star 

In this paper we present numerical results describing the structure of neutron star 
based on the relativistic mean field theory. It is possible to describe a static spherical 
star solving the OTV equation. 

dP(r) _ G P(r), (m(r) + jP(r)r 8 ) 



C 2 T 



dm(r 



4irr 2 p(r) (52) 
dr 

Having solved the OTV equation the pressure p(r), mass m(r) and density p(r) were 
obtained. To obtain the total radius R of the star the fulfillment of the condition 
p(R) = is necessary. This allows to determine the total gravitational mass of the 
star M{R). 

Introducing the dimensionless variable £, which is connected with the star radius r 
by the relation r = a£ enables to define the functions p(r), p(r) and m(r) 

p(r) = PoX(x(0) (53) 
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Figure 7: The neutron star density profile for fermions pp scalar bosons <p and gauge 
bosons (u^pfy. 

P(r) = P ip(x(t)) (54) 



m(r) = M Q v(£) (55) 

by £. If we define dimensionless functions. 

GM Q p c M 4 3 

= -p^-. Ai = 3^, M c = - W (56) 

are also need to achieve the OTV equation of the following form 

= " A(x(0 + y ®> ^(i _ asgl) (57) 

| = /*(0* 2 (58) 

with r 9 being the gravitational radius. The equations ( p^JSEjj ) are easy integrated nu- 
merically, For example, for the neutron star with the central density p c = 9 10 14 g/cm? 
the star profile in the mean field approach is presented on the Fig.0. It's interesting 
that fermions contribution to density is lower then a half of total density. The biggest 
contributions come from the nucleons, the gauge boson field and scalar boson field 
<p. On the surface of star we have non zero value of density (p ~ 0.25 < p ) an d 
nucleon mass lower than vacuum nucleon mass (mj? ~ 762 Me V < M). Inside star 
with p c = 9 x 10 14 g/cm 3 (maximal stable configuration for the TM1 parameters set) 
value of nucleon mass varies from ~ 300 MeV in centrum of star to ~ 762 MeV on 
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Figure 8: The effective neutron mass profile (mjv = MS) inside the star. 




Figure 9: The R-M diagram for the neutron star. 



the surface (Fig. m. Fig. 10. shows the stellar masses as a function of the central 



density. The parameters of the maximum mass configuration are: 

M max = 1.91 M , R = 12.84 km. (59) 

This fact is easy to notice on the mass-radius diagram (Fig. H). When temperature 
is different from zero, a star is in generar bigger and more massive (see Figs. || and 

m. 



Conclusion 

The structure of static neutron stars can be determined by solving the Tolman- 
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Figure 10: The neutron star radius dependence from the central density p c . 

Oppenheimer-Volkoff equations. The equation of state will strong influence on the 
neutron star properties. The aim of this work has been to present the equation of 
state in the RMF for the nuclear matter with temperature different from zero. The 
hot neutron star in the simple L2 parameters set was examined in our previous work 



Ej]. The main objective of our work was to study the influence of the temperature 
on the main parameters of a neutron star. Neutron matter at finite temperature is 
of increasing interest in relation to the problems of hot neutron stars. In order to 
achieve the proper form of the equation of state the relativistic mean field approach 
was involved. 
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Appendix: Nuclear matter properties 

The nuclear symmetric matter may be described by the phenomenological equation 



of state 14 



e(u) = Pa u(M + |H + \Au- + ^u°) (60) 

where u = p/po is a dimensionless density. For symmetric nuclear matter dimension- 
less po = 2.5 10 14 g cm~ 3 = 0.15 nucleons / fm 3 = 140 MeV/m~ 3 . The parameter 
of e (|60| ) are 

K + 2E Ffi 

B = {^)[\E Ffl -E,l (62) 
a — i 6 

A = E - ^E F>0 - B. (63) 

Epfi is the nonrelativistic nucleon Fermi energy at the saturation point kp,o (see 
Table §). F° r the TM1 parameter set (see Table |l|) we have 

a = 1.5485; A = -1.58.52 MeV; B = 107.689 MeV; 

In this approach the binding energy (see Fig.[TI] ) is 

Eq(u) = E Ffi J + \Au + (64) 
2 (a + 1) 

The pressure is 

d(e/QB), 2^ s 1 . 9 Bo , . 

The minimum of the binding energy determine the equilibrium Fermi momentum 
kpfi, density po and incompressibility factor 

K = 9p — \ P0 . (66) 



The empirical value of K is 210 ± 30 MeV ||TTJ. The incompressibility factor is equal 
to 

K(u) = V , P \ Po = 10E FO u% + 9(Au + Bau u ). (67) 
du 

The chemical potential defined for particle species i is given by 

( \ ^ ffiQ\ 

= (68) 
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E/Q-M ( MeV 




Figure 11: The binding energy for nucleon symmetric TM1 phase according to the RMF 
approach (the solid line) and (|64|) (the dot line), respectively. 



v/c 



P'Po 



Figure 12: The adiabatic sound speed as the function of the relative density u = p/po for 
the TM1 parameters set. 



For example, for nucleon symmetric saturation point p = 912.73 Me V for nucleons. 
The adiabatic sound speed 

,V.n dP K , . 

- = -r = — 69 

V de 9/i v ' 

is presented on the Figjl^. At the saturation point we have the sound speed v = 
0.184c. For TM1 parameters, for p above 5p the theory looses its sense. For 
other parameters sets when K is grater the range of the theory validity is smaller. 
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